Dynamics of the entanglement rate in the presence of decoherence 
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Entanglement plays an essential role in quantum in- 
formation theory, the sharing of entanglement between 
sender and receiver allows for quantum teleportation[l], 
quantum superdense coding [2] and the other applications 
to quantum information processing [3]. Creating entan- 
glement in a proper way is thus an important issue. 

In general, entanglement between two systems can be 
generated if they interact in a controlled way. However, 
for a practical experiment, the production of entangle- 
ment is very difficult due to the weak interaction between 
the systems. Thus how to improve efficiency of the pro- 
duction by using those interactions become a very rel- 
evant problem. Very recently, Diir et al. [4] consider a 
situation that one has a given non-local Hamiltonian and 
ask, what is the most efficient way of entangling parti- 
cles? Their answers are that (i) the initially entangled 
two particles can improve the efficiency of the produc- 
tion, and (ii) one can also improve the efficiency by using 
some ancillas. 

In this paper, we shed light on this issue again by tak- 
ing the decoherence effects into account. As you will 
see, the problem for mixed states is complicated, thus we 
choose two special sets of mixed state to study the prob- 
lem. This is the limitation of this paper. The results 
show that there is a competition between entangling and 
decohering in the entanglement production process, some 
initial states may work very well in the absence of deco- 
herence, but they do not provide the best way to produce 
entanglement in the presence of decoherence. 

To begin with, we recall the definition of entanglement 
rate T{t)[A] 



T(t) = 



dE(t) 
dt ' 



(1) 



where E(t) denotes an entanglement measure of a state 
p(t). In this paper, we pay our attention first to the case 
of two qubits, and then generalize the discussion to the 
case of d-level system with d > 2. We use the following 
notations throughout this paper: = 1,2,3,4) stand 



for bases of the two-qubit system, = 1 00) , |^) = 
|01),|£ 3 ) = |10), |e 4 ) = 111)- \m){m = 0,1) denotes 
the two states of one qubit, and pij = (£i|p(i)|£j) rep- 
resent the matrix elements of p(t) in the space spanned 
by |&)(* = 1, 4). With those notations, we choose 
15 independent variables to describe a general two-qubit 
state p(t), they consist of 3 independent diagonal ele- 
ments Pn, /522> P33 an d 6 complex off-diagonal elements 
Pij(i,j — 1, ...,4, and j > i) of matrix p(i)[5]. 

In order to calculate the entanglement rate, we have 
to express the entanglement measure of the state p(t) as 
a function of pij(i,j = 1, ■■•,4, j > i). For an entangled 
two-qubit system, we may choose the Wootters concur- 
rence as the entanglement measure 



E(t) = S(c(p)), 



(2) 



with 



£(c) = K o )> 



h(x) — —xlog2X — (1 — x)log2(l — x), 



c(p) 



c{0, Ai - A 2 - A3 - A4}, 



where the A^s are the square roots of the eigenvalues of 
the non-Hcrmitian matrix pp with p — (a y a y )p*(a y (8 
a y ) in decreasing order. In the space spanned by i = 
1,2,3,4}, p reads 




-P34 

P33 
(P23)* 

-(pisY 



-P24 

P23 

P22 
-(Pl2)' 



Pli 
~P13 
-Pl2 

Pll 



(3) 



The Wootters concurrence gives an explicit expression 
for the entanglement of formation, which quantifies the 
resources needed to create a given entangled state. Note 
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that the Wootters concurrence is a function of Pij(i,j = 
1, ...,4, j > i). Therefore, we can write 



r(t) 



4 



8E dpi 



■ ■ i ^ dpi 



at 



, {Pij # P44)- 



(4) 



Eq.(4) shows that given a particular entanglement mea- 
sure E(pij), we just have to determine dpij/dt . In or- 
der to do that, we need to find the time evolution of 
the state p(t) in the presence of an environment. Gen- 
erally speaking, interactions between a quantum system 
and its environment result in two kinds of irreversible ef- 
fects: dissipation and dephasing. The first effect is due 
to the energy exchange between the system and its envi- 
ronment, whereas the second one comes from the system- 
environment interaction that does not change the system 
energy. Both dissipation and dephasing lead to decoher- 
ence. In what follows, we first drive an expression for the 
time derivative of p(t) in terms of Kraus operators, this 
equation may be useful in the case of the Krause opera- 
tors are easily given. Then we adapt the other description 
of decoherence to put forward our discussion. Consider 
a quantum system of two qubits p interacting with an 
environment pe — ^2 u Pv\v){v\, after a finite time evolu- 
tion governed by unitary evolution operator U(t,0), the 
total density operator(the system plus the environment) 
pt{t) is given as 



Pt(t) = U{t,0)(pE®p)U*(t,0). 



(5) 



Taking a partial trace over environment variables we can 
get the density operator of the two-qubit system in the 
following form [6] 

p(t) = Tv E p t (t) = K ^ 0)p(0)^(t, 0), (6) 



where K^(t, 0) = (p\^/p^U(t,0)\v}. The Kraus oper- 
ators K^ v satisfy ^ K^K^ = 1. No environment 
around the two-qubit system indicates that there is only 
one term in the sum eq.(6). In weak system-environment 
interaction limit, the density operator of the environment 
Pe remains unchanged in the whole time evolution pro- 
cess, this approximation can be found in the derivation 
of the master equation, which we will discuss later on. 
Under weak system-environment interaction, the Kraus 
operators can be expanded to first order of dt as 



K^(t + dt,t) = <ji\y/j£U(t + dt,t)\u) 

— s/pZSfiv - idt^/p^{p\H t \v) 



(7) 



where H t stands for the Hamiltonian of the total 
system(two-qubit plus its environment), and this equa- 
tion holds only for very small dt. Substituting Eq.(7) 
into Eq.(6), we obtain p(t + dt) in terms of p(t)(to first 
order of dt) 



pit + dt) = P (t) - idtj2v^(v\ H tW)p(t) 

fl.V 

+ idt P (t)j2vp^(^\HtW) 

= p(t) + dp(t). (8) 

In order to calculate the entanglement rate, we have to 
compute dpij/dt. Using standard perturbation theory, 
we find dpij/dt as follows 



with 



dpi 3 
dt 



d P (t) 
dt 



= (^\dp(t)/dt\^) 



-i{j2vp^(v\ H tW),p(t)} 



(9) 



Eq.(8) and Eq.(9) show that v ^/pZ{p\H t \v) play a 



role of effective Hamiltonian for the two-qubit system, in 
this sense we may rewrite Eq.(8) in the following form 



(10) 



where H e = £^ v ^/p^{p\H t \v) . This expression is useful 
and easy to handle when we know the Kraus operators. 

The other tool to study the quantum dissipative sys- 
tem is the master equation, which can be obtained in 
Markovian limit [7-9] . This approximation is very useful 
because it is valid for many physical relevant situations 
and its numerical solutions can be easily found. As given 
by Gardiner, Walls and Millburn, Louisell in their text- 
book [9], the reduced density matrix p of the open sys- 
tem which is linearly coupled to its environment obeys 
the following master equation of Lindblad form [10] 

p(t) = -i[H ,p] 



+ lYl G m (2X+pX- X-X+p pX~X+) (11) 



with 

K m = 2Rc 



G v 



2Re 



dTe^ T Ti env {A m (r)Al (0)p env } 



(-00 

/ dTe^ T Ti env {Al(r)A m (0)p env } 
Jo 



Here, p(t) = p(t, K m ,G m ) stands for the density opera- 
tor of the system and p env denotes the density operator 
of the environment, X^ are eigenoperators of the sys- 
tem satisfying [Hq,X^\ = ±huj m X^ l , H stands for the 
free Hamiltonian of the system, and A m {A\ n ) are oper- 
ators of the environment through which the system and 
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its environment couples together. Notice from Eq.(ll) 
that G m should vanish at zero temperature T = 0, while 
K m should not if A m are indeed destruction operators of 
some kind. 

The time derivative of the matrix elements pij in this 
case is 



dt 



dt 



mm- 



(12) 



Further more, we consider a case of a two-qubit sys- 
tem coupling to environments that consists of a set of 
harmonic oscillators. In this case, the master equation 
takes the following form at zero temperature 

P (t) = -i[H s p-pH s ] 

2 

+ 1 ^(2aLp4 - <°-P - P°>-), ( 13 ) 
i=i 

where H s is the system Hamiltonian, which governs time 
evolution of the two-qubit system in the absence of its en- 
vironment, o l + {(j l _) are pauli matrices, and 7 represents 
the damping rate. Considering a system Hamiltonian 

H a = —^ + —^ + Kg{p\<ri+<jl<r l + ) (14) 
and substituting Eq.(13) into Eq.(12), we obtain(setting 

n = 1) 



Pll = 7(P22 + P3d), 

P22 = -igP32 + igp23 + 7P44 - 7P22, 
P33 = ~igP23 + igP32 + 7P44 - 7P33, 
P44 = -27P44, 

P12 = igpi3 + 7P34 - O.57P12 + iwpi2, 
P13 = igPu + 7P24 - O.57P13 + iwp 13 , 

P\A = -7P14 + 2iwpi4, 
P23 = ~igP33 + igP22 - 7P23, 
P24 = -igP3A - 1-57/924 + ^P24, 
P34 = -igP2A - 1-57/034 + iu>p3A 



is a family of such initial states. In terms of Bell bases, 
eq.(16) can be written as 

po = a\^+){^+\ + b\*i>-)(i>-\ + c|0+>^+i + d\r){r\- 

This is just the Werner state and any two-qubit entan- 
gled state can be expressed in this form by performing 
a random bilateral rotation on each shared pair [3]. We 
may use this kind of entangled state to demonstrate how 
it decoheres, although it could not be used to increase en- 
tanglement. We would like to mention that for a general 
mixed state characterized by 15 independent parameters, 
the problem becomes complicated. Hence we choose two 
special sets of initial states to get some insights into the 
formalism. In terms of Wootters concurrence, the entan- 
glement measure for states(16) is 

E(t) = -xlog 2 x - (1 - x)log 2 (l - x), (17) 

where x = (1 + \/I - F 2 )/2, F = a - b - c - d for 
a — b — c— d > 0, or for a — 6 — c — d < 0, and we 
assume a > b,a > c,a > d without loss of generality. 

It is easy to show that the entanglement rate defined 
by (1) takes the following form(a + b + c + d = 1) 

dE(t) da 



T(t) 



da dt ' 



(18) 



with 



dE(t) 
da 
da 

dt 



= (log 2 x - log 2 (l - x)) 



F 



VI -F 2 ' 



c + d 



-7 — aj. 



(19) 



(15) 



By definitions, F > and x > |, then §f > 0. So, 
wether the entanglement increase or decrease for the ini- 
tial states (16) only depends on And ^| is always 
below zero. Therefore, we could not increase entangle- 
ment starting from the initial states (16). Because ^ 
does not depend on c and d, the entanglement rate de- 
pends on c and d linearly, and with c + d increase, the 
entanglement decreases. Whereas the entanglement rate 
is inversely proportional to 7. The dependence of T(t) 
on parameter a is shown in figure 1. 

Rate 



The Hamiltonian _ff s (14) describes two two-level atoms 
with dipole-dipole interactions, which are a source of 
creating entanglement for trapped atoms in an optical 
lattice[ll,12]. For an initial state, if the interaction 
terms(with coupling constant g in Eq.(15)) have no ef- 
fects in the time evolution process, it (an example is given 
below) could not be used to create entangled state or to 
increase entanglement. Initial state with the form of 
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(16) 
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FIG. 1. Dependence of the entanglement rate T on a with 
7 = 0.01 and c + d = 0.1. 
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This figure shows that the larger the parameter a, the 
smaller the entanglement rate. For a limit case a = 1, 
and c = d = b = that corresponds to the maximally en- 
tangled state the entanglement rate takes its min- 
imum over states (16). 

In contrast to example eq.(16), we present here another 
kind of states 



Po = 



'0 0> 

p q 

q* 1-p 

,0 0, 



(20) 



This kind of state is of interest because entanglement con- 
tained in those states range from zero to one(maximally 
entangled state) . And it is a typical family of states for 
a ring of N qubits in a translation invariant quantum 
state[13]. 

The entanglement measure for this family of states has 
the same expression as Eq.(17) except for replacing F by 
G 

G = 2\q\. (21) 

The positivity of the state po Eq.(20) require 

R = p 2 -p+\q\ 2 <0, (22) 

this indicates that \q\ < \ for all family of states (20). 
And for \q\ = i, there is only one value available for 
p, i.e. p = This is shown in Fig. 2 which gives the 
dependence of R on p and \q\. 



dE(t) _ 
dp 

dp r 1 , 

f t = -2gq J + - 7 (1 - 2p), 

8E(t) G dG 

—^j- = \ l °92V - log2{l - y))- 



VI - G 2 dq 1 ' 



dg 1 
~dt 
dE(t) 



= g(2p - 1) - jq 1 , 
= {log 2 y - log 2 {l - y)) 



G dG 



VI - G 2 dq 1 



8q R R 
dG 4q R W 



dqW) \q\ 



(24) 



where y = |(1 + V1 - G 2 ). Equations (23) (24) together 



give 



r(t) = io 92 



y 



G Agq I (2p~l)^A 1 \q\ 



l-y VT - G 2 ' 



This shows that there are competitions between entan- 
gling and decohering. If g/j > \q\ 2 /(q I (2p-l)), T(t) > 0, 
entanglement increases. Otherwise the entanglement de- 
creases. In other words, in order to get a positive entan- 
glement rate, the decoherence rate 7 and the coupling 
constant g should satisfy condition g > j\q\ 2 / (q 1 (2p — 1)) 
For parameters p = 0.6, g = 0.2, 7 = 0.01, the entangle- 
ment rate F versus q 1 and q R is illustrated in Fig. 3. 



0.2 




0.5 L 



FIG. 2. Dependence of the quantity R on p and \q\, that 
shows which region of |g| and p axe available for the state(20). 

The entanglement rate in this example is 

dE(t) dp dE(t) dq R dE{t) dq 1 



T(t) 



dp dt 



l R dt dq 1 dt 



(23) 



where q R ^ represents the real (imaginary) part of q. 
Eq.(15) and (17) together give 



Rate 




FIG. 3. Entanglement rate for the family of state 
(20) versus q 1 and q R . The parameters chosen are 
p = 0.6, g = 0.2,7 = 0.01. 



The maximum of the entanglement rate is 0.4 corre- 
sponding q 1 — 0.5, q R = 0.0, while the minimum of T is 
-0.1 at about q 1 = 0, q R = 0.5. Similarly, the maximum 
and minimum of the entanglement rate change with g 
and 7, but ^Lztmm)' 9ma:r( m zn) which corresponding to 
- 1 - max (Tmin) do not. 
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In the end of this paper, we generalize the formulas de- 
rived above to the case of multilevel systems, we denote 
by Ui the generators of the group SU(N) with N being 
the dimension of the Hilbert space of the system A, and 
Ti the generators corresponding to system B with dimen- 
sion M. With this notation, we may write any density 
matrix of the composite system A and B in the following 
general form 



(25) 



Here we choose on, (3i and jij as the independent vari- 
ables to characterize the state of system A plus B. The 
entanglement measure E(on, /3j,7»j) is thus a function of 
a.i , (3j , and 7^ . It is natural to express the entanglement 
rate in the following way 



r(t) 



EdE don y—^ dE dfli y—^ 



dE dpjj 
djij dt 



(26) 



Given an entanglement measure E(cti, (3j, 7^), the 
derivatives dE/don, dE/dfij and dE/djij can be easily 
calculated. The remained task is only to determine the 
time derivative of at, flj, and 7^. Proceeding as before, 
we obtain 



d/3i dpAsif) 
— =Tr B (r i —^), 



dt 



(27) 



Where dpAs/dt has a similar expression with Eq.(lO) 
or Eq.(ll), depending on what formalism you choose to 
describe the time evolution of the system. 

In summery, taking the dccohercncc effects into ac- 
count, we study dynamics of the entanglement rate for 
two special sets of initial states. The interaction un- 
der consideration is of pairwise. The results show that 
there are competitions between decohering and entan- 
gling, those competitions lead to (l).For a specific inter- 
action and a dccohercncc scheme, some initial state could 



not be used to prepare entanglement. (2). Some initial 
states can be used to prepare or increase entanglement 
under a proper choice of the parameters. 
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